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Abstract: We have investigated the resonant production of a stop at the Large Hadron 
Collider, driven by baryon number violating interactions in super symmetry. We work in the 
framework of minimal supergravity models with the lightest neutralino being the lightest 
supersymmetric particle which decays within the detector. We look at various dilepton 
and trilepton final states, with or without b-tags. A detailed background simulation is 
performed, and all possible decay modes of the lighter stop are taken into account. We 
find that higher stop masses are sometimes easier to probe, through the decay of the stop 
into the third or fourth neutralino and their subsequent cascades. We also comment on 
the detectability of such signals during the 7 TeV run, where, as expected, only relatively 
light stops can be probed. Our conclusion is that the resonant process may be probed, at 
both 10 and 14 TeV, with the R-parity violating coupling Ag 12 as l° w as 0.05, for a stop 
mass of about 1 TeV. The possibility of distinguishing between resonant stop production 
and pair-production is also discussed. 
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1. Introduction 

The current structure of the standard model (SM), with gauge invariance and renormal- 
isability built in, implies automatic lepton and baryon number conservation. This is no 
longer true in the supersymmetric (SUSY) extension of the SM 111 |2| , where scalars carry- 
ing baryon or lepton number are present. Thus the superpotential of the minimal SUSY 
standard model (MSSM), namely 



Wmssai = hfjQiD'Hd + h^QiU^H u + //!,/.,/•:; 7/,, + fiH u H d (1.1) 
can in principle be augmented to include 



y^RPV = HiLiH u + \ij k LiLjEl + \' ijk LiQjDl + X'^UfD^Dl (1.2) 

which contain terms that are gauge invariant and renormalisable but explicitly violate 
lepton or baryon number. Here, L(E) is an SU(2) doublet (singlet) lepton superfield and 
Q (U,D) is (are) an SU(2) doublet (singlet) quark superfield(s). H u and are the two 
Higgs doublet superfields, /i is the Higgsino parameter and (i, j, k) are flavour indices. Each 
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term in equation (1.2) violates R-parity, denned as R = (— l) 3 ( B ~ i )~ 2S ' (where B is baryon 
number, L is lepton number and S is spin), against which all SM particles are even whereas 
all superpartners are odd. The consequence of violating R-parity is that superpartners 
need not be produced in pairs anymore, and that the lightest superparticle (LSP) can now 
decay. The strongest argument for studying R-parity violation is that it does not arise as 
an essential symmetry of MSSM. However, the requirement of suppressing proton decay 
prompts one to allow only one of B and L to be violated at a time. 

The collider phenomenology in the absence of R-parity may be very different from that 
of the usual R-parity conserving MSSM. In particular, if the R-parity violating (RPV) 
couplings are large enough, the LSP will decay within the detector and one no longer 
has missing- Ifrp sis ci convenient discriminator. Although studies have taken place on such 
signals, closer looks at them are often quite relevant in the wake of the Large Hadron 
Collider (LHC). In particular, it is crucial to know the consequences of broken R-parity 
in the production of sparticles. Here we perform a detailed simulation in the context of 
the LHC, highlighting one possible consequence of the B- violating term(s), namely, the 
resonant production of a squark — in this case, the stop. 

Many of the RPV couplings have been indirectly constrained from various decay pro- 
cesses, including rare and flavour- violating decays and violation of weak universality. The 
constraints derived are of two general kinds — those on individual RPV couplings, assuming 
the existence of a single RPV term; and those on the products of couplings when at least 
two terms are present, which contribute to some (usually rare) process. The constraints 
obtained so far are well-listed in the literature (3[. 

The L- violating terms are relatively well-studied, partly because of their potential role 
in generating neutrino masses and are constrained by indirect limits. In comparison, the 
baryon-number violating coupling are relatively unconstrained. X'( 12 113 are constrained 
from double nucleon decay and neutron- antineutron oscillations |!| The rest of the 
couplings are constrained only by the requirement that they remain perturbative till the 
GUT scale. Limits on A^ • type of couplings due to the ratio of Z-boson decay widths for 
hadronic versus leptonic final states have been calculated for a stop mass of 100 GeV |7j. 
However, the results do not restrict the couplings for high stop masses of concern here. 
The coupling A3 is thus practically unconstrained for large stop masses. It is also known 
that mixing in the quark sector causes generation of couplings of different flavour structures 
and can therefore be constrained by data from flavour changing neutral currents(FCNC)||. 
Such effects arising from mixing in the quark and squark sector can affect the contribution 
of R-parity violation to physical process and alter the limits |J. However these effects are 
model dependent and have not been taken into account here. 

It has been already noticed that such large values of A"-type couplings as are still al- 
lowed, not only cause the LSP to decay, but also lead to resonant production of squarks via 
quark fusion at the LHC. The rate of such fusion can in fact far exceed that of the canon- 
ically studied squark-pair production. One would therefore like to know how detectable 
the resonant process is at the LHC. Furthermore, one needs to know the search limits in 
different phases of the LHC, and how best to handle the backgrounds, both from the SM 
and the R-conserving SUSY processes. These are some of the questions addressed in this 
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paper. 

Single stop production, mostly in the context of the Tevatron, was studied in detail 



in[lC, 11]. A full one-loop production cross section can be found in[12|. A study of SUSY 
with the LSP decaying through baryon-number violating couplings and therefore giving 
no missing energy was done in |R|]. Further studies on determining the flavour structure 
of baryon number violating couplings and possible mass reconstruction following specific 
decay chains can be found in|l4|, |15[|. There have also been recent studies on possible 



LSPs|l6[ and identification of R-parity violating decays of the LSP using jet substructure 



methods |17|. A recent study on identification of stop-pair production via top-tagging using 



jet-substructure can be found in [18|. 



We find, however, that the earlier studies on resonant stop production are inadequate 
in the context of the LHC. We improve upon them in the following respects: 

• In the work done for the Tevatron, the sparticle masses were required to be less 
than 500 GeV to be within reach. Thus, the gluino was also required to be much 
lighter than a TeV to avoid large radiative corrections to squark masses. This implied 
that in the constrained MSSM (cMSSM) [19] scenario, the LSP, assumed to be the 
lightest neutralino (x5)> had to have mass less than 100 GeV. If we allow only the 
term proportional to A^,-, the only three-body decay of \\ is X? ~~ ^ ididj(tds). Since 
the neutralino is much lighter than the top, it can decay only via a 4-body decay and 
therefore has long lifetime and decays outside the detector for all allowed values of 
Aa,„ [pOj. Thus, one still has the canonical missing- Et signature. This was one of the 



main assumptions in [ 11 1 . However, if the stop mass is beyond the Tevatron reach 
but within the reach of the LHC, we may indeed have lightest-neutralino mass high 
enough to allow decay within the detector. 

We focus on this richer and more challenging scenario within the framework of min- 
imal supergravity (mSUGRA) models. We investigate the LHC reach for detection 
of the lightest stop assuming that the lightest neutralino is the LSP which decays 
within the detector and the stop is heavy enough to be beyond the reach of Tevatron. 

For a light stop, the only available R-parity conserving decay modes are into the 
lightest neutralino (x?) and lighter chargino(x^) i-e. t — > txi,bxf- For stop mass 
near a TeV, the decay modes into higher neutralinos and the heavier chargino may 
open up, leading to different final states. We have found that this drastically improves 
the detectability of the signature over the SM backgrounds. 

We have taken into account all the potential backgrounds at the LHC, including those 
from tt + jets, Wti + jets, Zti + jets, which pose little problem at the Tevatron. A 
detailed investigation towards reducing these backgrounds has been reported in the 
present study. 

In identifying signals of stop decay, it is often helpful to tag b quarks. This is 
especially important since the reconstruction of energetic top quarks in multi-top 
final states is difficult and has a low efficiency. However, the b-jets produced from 
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stop decay can be quite hard, especially when the stop is heavy and one is looking 
at the decay in the bxf channel. It is not clear that the b-tagging efficiency is 
appreciable at the LHC for b-jets with p? 100 GeV. With this is mind, we have 
performed a conservative analysis with a b-jet has zero detection efficiency unless its 
Pt lies in the range 50 - 100 GeVpI|. 

We present our results for centre of mass energies of 7, 10 and 14 TeV at the LHC. 
In section 2, we discuss the rates for resonant stop production, the different decays of 
the stop and our choice of benchmark points to account for all of them. In section 3, we 
present a detailed description of cuts required to isolate the signal and section 4 contains 
the numerical results from our simulation. We also comment on the distinguishability of 
such a signal from dilepton signals coming from R-parity conserving MSSM. It is possible 
to have LSPs other than Xi when R-parity is violated. We comment on the possibility of 
detection in such cases in section 5 and summarise and conclude in section 6. 



2. Resonant stop production and decays 
2.1 Stop production 

The resonant stop production process depends on B-violating couplings proportional to 
Agj -, and also on fraction of the right-chiral eigenstate in the mass eigenstate concerned. 
We concentrate on the production of i\ since the lighter stop eigenstate usually has a higher 
fraction of t^. The resonant production cross section is given by 



2vr sin 2 9j 
fl 3m? 

dxxte-Miixdfjix*) + fifaWxWQ. - -±) (2.1) 

where sin^ is the amplitude of finding a tn in t\, fi is the proton parton distribution 
function for a parton of species i and £(1,2) are the momentum fractions carried by the 
respective partons. Out of the three possible A" couplings, contributions via and A323 
are suppressed due to the small fraction of b quarks in the proton. We therefore look at 
the production of top (anti) squark through the fusion of the d and s (anti)quarks, via the 
coupling A 3 ' 12 . Since the actual cross section for production of the lightest stop depends on 
the mixing angle via sin 2 9i, it is useful to define the cross section in terms of an effective 
coupling Agjj = sin 8 '^-^312 ■ 
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The production cross section at the LHC with centre-of-mass energies of 7, 10 and 
14 TeV is given in Figure |l|. As an illustration, we have chosen the value X"^ = 0.2 
which is consistent with the existing limit on A^- In general, both t\ and ti will be 
produced. However, due to larger mass and smaller fraction of tn , t<i is rarely produced. For 
comparison, we also present the ii-pair production cross-section via strong interaction. For 
> 500 GeV, the resonant production dominates over pair-production for > 0.01 



rriT 



at 14 TeV. Resonant production can therefore hold the key to heavy stop signals if baryon 
number is violated. 



14 TeV 
10 TeV 
7 TeV 
14 TeV 
TX 10 TeV 
\\ 7 TeV 




400 



800 1200 
Mr (GeV) 



1600 



2000 



Figure 1: Production cross section at the LHC for <fs = 7, 10 and 14 TeV with \'L f = 0.2. The 
corresponding cross sections for R-conserving pair production are also shown. 



At next-to-leading order, the production cross section at yfs = 14 TeV is modified 



by a k-factor of about 1.4[12|. The uncertainty due to renormalisation and factorisation 
scales at lowest order is about 10% and drops to 5% when NLO corrections are taken into 
account. 



2.2 Stop decays and choice of benchmark points 

We wish to make our conclusions apply broadly to a general SUSY scenario and to include 
all possible final states arising from stop decay. However, the multitude of free parameters 
in the MSSM often encourages one to look for some organising principle. A common 
practice in this regard is to embed SUSY in high-scale breaking scheme. Following this 



practice, we have based our calculation on the minimal supergravity (mSUGRA) model] 22]. 
mainly for illustrating our claims in a less cumbersome manner. The high scale parameters 
in this model are: mo, the unified scalar parameter, mi/2, the unified gaugino parameter, 
sign(fi), where fi is the Higgsino mass parameter, Aq, the unified trilinear coupling and 
tan /?, the ratio of the two Higgs vacuum expectation values. 

Although the production cross section of the stop depends only on the mass and mixing 
angle of the stop, any strategy developed for seeing the ensuing signals has to take note 
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Figure 2: Lighter stop decay branching fractions in different modes for tan/3 = 5, Aq = —1500 
(top left) ; tan/3 = 40, A Q = -1500 (top right) and tan/3 = 10, A = (bottom). A^' 12 = 0.2, \i > 
and m 1 / 2 — 450 GeV in all cases. 



of the decay channels. We have tried to make our analysis comprehensive by including 
all possible decay chains of the stop. Thus we have included decays into tx® , bxf, tg and 
ds, of whom the first three are R-conserving decays while the last one is R-violating. The 
charginos, neutralinos or the gluino produced out of stop-decay have their usual cascades 
until the LSP (here Xii the lightest neutralino) is reached. The X\ thereafter undergoes 
three-body RPV decays driven by X'^ 12 , to give rise to final states consisting leptons and 
jets of various multiplicities. 

We observe that for the same values of (mcm^), the mass and branching fractions 
of the stop may vary drastically with different values of (tan /3, Aq). We shall choose /i > 
for all the benchmark points as it is favoured by the constraint from the muon anomalous 
magnetic moment |2^j . 

Since we explicitly want to study the situation in which the neutralino decays within 
the detector, the only available decay mode is Xi ~^ tds(ids). We therefore require that 
the neutralino mass be greater than the top mass to allow for a three-body decay. We 
fix mi/2 = 450 GeV which gives M^o ~ 180 GeV. We also choose the high scale value of 
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A312 ~ 0.065 such that it gives a value of 0.2 at the electroweak scale. 

Figure ^ shows the branching fractions into various final states for three different 
choices of (tan j3, Aq), namely, (5, —1500), (40, —1500) and (10, 0) for different stop masses, 
obtained by varying uiq. We notice that, for low mo, the dominant decay mode is bx\ in 
the third case of Figure [2], while it is t~x\ in the first two cases. We also notice that the 
decays into higher neutralinos and charginos open up earlier for tan f3 = 40 and compared 
to tan j3 = 5. 

The Tevatron reach for single stop production is about 450 GeV. We therefore start 
with a benchmark point with stop mass of 500 GeV, just beyond this reach (Point A). The 
major decay channels in this case are txi,bxf- A stop mass of a TeV at the electroweak 
scale may be obtained by various configurations in the high-scale parameter space. How- 
ever, from the above plots, one expects its decays to change significantly with different 
parameters. Our objective is to determine whether signal of resonant production of a stop 
of mass near a TeV can be probed irrespective of what the high-scale parameters are. For 
this, we fix Mi ~ 1 TeV. We first look at the case with Aq = —1500. We construct 
two benchmark points with tan j3 = 5(Point B) and 40(Point C) which correspond to the 
opposite ends of the allowed range in tan /3. We see that for a stop mass of 1 TeV, the 
decays into the higgsino-like xt an< ^ X3 become dominant goes to high tan /3. 

Similarly, we also look at a point with Aq = tan/3 = 10 (Point D). In this case, 
we find that the Higgsino channels open up fairly early and the dominant decay is 
followed by t\%. As we shall see in the next section, this plays a crucial role in enhancing 
multi-lepton signals of a resonantly produced stop. Finally, since the decay into a top and 
a gluino does not open up until much higher stop masses, we also construct one point in 
which the stop decays dominantly into tg (Point E). 

Points A, B and E correspond to the same value of (tan/3, .Ao) = (5,-1500) and 
therefore provide a description of how the signal changes when only thq is varied. This 
choice of parameters also corresponds to the most conservative case in terms of signal since 
the decay modes into the higher gauginos does not open for a large region in the parameter 
space. We will therefore use these points to obtain limits on 

We have tabulated the parameters and significant decay modes in Table [l|. The bench- 
mark points were generated with RPV renormalisation group running of couplings and 
masses using SOFTSUSY 3.0.2[§4| and the RPV decays were calculated with the ISAWIG 
interface to Isajet[p5]1 . 

The decay width of the stop in the R-parity violating channel ds depends only on 
Ag^-j and the stop mass. Therefore, the branching ratio into this channel for same values 
of Agjj and stop mass depends only on the decay widths of the other channels open at 
the same time. For the benchmarks under consideration, Xi 2 an d Xi nave large gaugino 
fractions whereas X3 4 and X2 have large higgsino fractions. The large top mass means 
that stop coupling to higgsino-like chargino and neutralinos is large. Thus as soon as these 
decays become kinematically allowed, they quickly dominate over the decays into gaugino- 
like chargino and neutralinos. This can be seen for points B, C and D which have nearly 
identical stop masses and sin 2 Q~ t (Ag 12 = 0.2 at electroweak scale for all points). Large 
tan/3 opens up the bx£ mode early in point C as compared to point B and makes the 
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(m ,tan(3,A ) 


Mr 
*1 


sin 2 6~ t 


Dominant decay modes 
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t\i (U.60), OXi (U.iyj, CIS (U.48J 


B 


(1650,5,-1500) 


1002 


0.97 


txi (0A8);tx° 2 (0.04); bxf (0-10); ds (0.38) 


P 


( 1 ^7n /in 1 ^nn"i 


1 nno 

1UUZ 


u.yo 


IXi (,U.dDJ,tX 2 (.U.UDJ, 

bxf (0.12); 6x? (0-21); ds (0.27) 


D 


(1250,10,0) 


1008 


0.97 


txS (0.13); txl (0-04); 4*° (0-20); t X ° (0.13); 
bxt (0.08); 6xJ (0.33); ds (0.10) 


E 


(2450,5,-1500) 


1404 


0.99 


^ (0.39); ix? (0.15); t X ° 2 (0.02); **§ (0.08); 
txi (0-07); b*+ (0.02); 6 X + (0.17); ds (0.11) 



Table 1: Benchmark points and the dominant decay modes of the lighter stop. = 0.2, /i > 
and mi/2 — 450 GeV for all benchmark points. 



branching fraction into ds for point C much lower. For point D, the branching fraction 
into higgsino-like chargino and neutralinos is larger than 60% and the RPV decay fraction 
is only about 10%. The t — t — g coupling comes from strong interactions and therefore the 
tg channel dominates whenever it becomes kinematically allowed (as in point E). 



3. Event generation and selection 
3.1 Event generation 



Signal events have been generated using HERWIG 6.510[26], and jets have been formed 



using anti-A;T algorithm[27] from FastJet 2.4.1. SM backgrounds have been calculated 



using Alpgen 2.13 p8|j showered through Pythia [2S] with MLM matching. We have used 



CTEQ6L1 parton distribution functions |3(|. The renormalisation and factorisation scales 
have been set at the lighter stop mass (Mj ) for signal, while the default option in ALPGEN 
has been used for the backgrounds. 

In R-parity conserving MSSM, the production of two heavy superparticles requires a 
large centre-of-mass energy at the parton level. This allows us to further suppress the SM 
background by applying cuts on global variables like the "effective mass" (M e ff). Since we 
no longer have a large missing- Er and the energy scale of the resonant production process is 
not very high, the SM background cannot be suppressed so easily. We therefore concentrate 
on leptonic signals with or without b-tags to identify the signal over the background. 



3.2 Event selection 

Decay of the lighter stop in this scenario can lead to a variety of final states. Out of them, 
we have chosen the following ones: 

• Same-sign dileptons: SSD 

• Same-sign dileptons with one b-tagged jet: SSD + b 

• Trileptons: 3/ 
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We do not consider the RPV dijet channel as a viable signature due to the enormous 
background from QCD processes. Similarly, we also omit opposite-sign dileptons due to 
large backgrounds from Drell-Yan, ti etc. 

We have imposed the following identification requirements on leptons and jets: 

• Leptons: A lepton (I) is considered isolated if (a) It is well separated from each jet 
(J): ARij > 0.4, (b) The total hadronic deposit within AR < 0.35 is less than 10 GeV. 
We consider only those leptons which fall within \rj\ < 2.5 with pr > 10 GeV. Here, 
AR = \J At] 2 + A<fi 2 where r\ is the pseudo-rapidity and <fi is the azimuthal angle. 

• Jets: Jets have been formed using the anti-for algorithm with parameter R = 0.7. 
We only retain jets with px > 20 GeV and \r/\ < 2.5. 

• B-tagged jets: A jet is b-tagged with probability of 0.5 if a b-hadron with 50 < pt < 
100 GeV lies within a cone of 0.7 from the jet axis. We have set the identification 
efficiency to be zero outside this window, in order to make our estimates conservative. 

We also apply the following extra cuts on various final states to enhance the signal over 
background: 

• Cut 1: Lepton-py: We demand that the pt of the leptons be greater than (40, 30) 
GeV for dilepton and (30, 30, 20) for trilepton channels. This cut removes the back- 
ground from semileptonic decays of b quarks. It strongly suppresses the bb + jets, 
Wbb + jets and ti+jets background in the SS'D-channel coming from semileptonic 
b-decays. 

• Cut 2: Missing Et- At least one lepton in the signal always comes from the decay 
of a W boson and is accompanied by a neutrino. We demand a missing- Et greater 
than 30 GeV from all events. This helps in reducing the probability of jets faking 
leptons. Missing- Et has been defined as \pr, V isible\- 

• Cut 3: Jet pt- We demand that the number of jets, rij > 2 with ptUi) > 100 GeV 
and pt(j2) > 50.0 GeV for SSD and SSD + b. This cut is useful when high stop 
mass is very high and the production cross section is very low. 

• Cut 4: Dilepton invariant mass: We also apply a cut on dilepton invariant 
mass {M1112) around Z-mass window (\Mi lt i 2 — Mz\ < 15.0 GeV) for opposite sign 
dileptons of same flavour in trilepton events. This serves to suppress contribution 
from Zbb + jets and WZ + jets background to trileptons. 

Due to the Majorana nature of neutralinos, A"-type interactions result in equal rates 
for tds and ids-type final states. Therefore, the most promising signals are those involving 
same-sign dileptons (SSD). This not only applies to Xi but also to the higher neutralinos 
produced in stop decay, whose cascades can give rise to Ws. SSD have previously been used 



extensively for studying signals of supersymmetry [31, ^] . The most copious backgrounds 



to SSD processes come from the processes tt and Wbb due to one lepton from W and another 
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from semileptonic decays of the b-quark. There is also a potentially large contribution from 
bb due to B° — B° oscillations along with semileptonic decays of both B-mesons. The effect 
of oscillations is simulated in the Pythia program. The pr-cuts on leptons have been 
selected to minimise the background from heavy flavour decays |33| . We find that after the 
isolation and pt cuts on leptons, Wbb and bb cross sections fall to sub-femtobarn levels. 

We simulate the ti + jets background up to two jets. The trilepton channel has 
another source of backgrounds in W Z + jets; however, we have checked and found them to 
be negligibly small after applying all the cuts. We also generate Wti + jets and Ztt + jets 
up to one jet. 

It should be mentioned here that the dilepton and trilepton final states can also arise in 
the same scenario from the pair-production of superparticles. These include, for example, 
pair production of gluinos and electroweak production of chargino-neutralino pairs. Such 
contributions have been explicitly shown in the plots in section 4. 

We also expect that the pt distribution of the t\ becomes significantly harder if the 
NLO corrections are taken into account [12]. Our cuts on leptons have been designed to 
cut off the background from semileptonic b-decays by requiring the pt to be about half 
the mass of the W. Therefore, if only the lepton cuts are used, we do not expect a large 
change in the efficiency of the cuts quoted in the next section. 



4. Results 

We present results using X'^ 12 = 0.2; the predictions for other values of this coupling can be 
obtained through scaling arguments. If is scaled by a factor of n then the production 
cross section as well as the decay width of the RPV channel scale by n 2 . All other decay 
widths remain unchanged. If / is the branching fraction of t — )■ ds before scaling, then the 
signal rates in any other channel are scaled by a factor of 

Rnew ^ /, ,\ 

iUT~ (n2-l)/ + l [ ■ > 

4.1 Limits at y/s =14, 10 TeV 

The numerical results for various signals corresponding to the five benchmark points for 
LHC running energy \fs = 14, 10 TeV are presented in Tables ||(for the SSD channel), || 
(for SSD + b) and | (for 3/). 

We can make the following observations from the numerical results: 

• The final states SSD and SSD + b consistently have substantial event rates at both 
14 and 10 TeV. Furthermore, the simultaneous observation of excesses in the SSD and 
SSD +b channel can serve as definite pointer to the production of a third generation 
squark. 

• For point A, which is just above the Tevatron reach, we can achieve more than 5a 
significance in the SSD channel with just 100 pb _1 data at both 14 and 10 TeV. For 
point E, which has = 1500 GeV, we can reach 3<r with 1(3) fb~ l and ha with 3(9) 
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26.3 


24.7 


10.0 


307.0 
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7.0 


3.6 


Wti + nj 


17.0 


9.2 


8.7 


5.2 


7.6 


3.9 


3.7 


2.0 


Zti + nj 


12.7 


6.7 


6.7 


4.1 


4.9 


2.3 


2.2 


1.4 


Total 


717.6 


42.2 


40.1 


19.3 


319.5 


14.9 


12.9 


7.0 



Table 2: Effect of cuts on signal and SM background cross sections (in fb) in the SSD channel at 
y/s = 14, 10 TeV. Cut refers to all events passing the identification cuts. Cuts 1-3 are described 
in the text. The numbers corresponding to best significance (s/^/b) of the signal (s) with respect 
to the background (b) are highlighted in bold. 



SSD + b 




14 TeV 








10 TeV 






Point 


Cut 


Cut 1 


Cut 2 


Cut 3 


Cut 


Cut 1 


Cut 2 


Cut 3 


A 


243.2 


134.7 


121.2 


14.1 


158.7 


67.5 


61.8 


6.6 


B 


13.8 


9.6 


9.2 


4.7 


8.3 


6.0 


5.6 


2.8 


C 


25.3 


15.2 


14.6 


8.0 


12.7 


7.8 


7.4 


4.1 


D 


47.9 


23.9 


23.0 


15.2 


21.2 


9.9 


9.5 


6.3 


E 


11.3 


6.2 


6.0 


5.3 


4.1 


2.2 


2.1 


1.8 


ti + nj 


173.0 


7.6 


7.1 


2.3 


80.9 


4.2 


1.4 


1.3 


Wti+nj 


6.7 


0.8 


0.8 


0.6 


4.0 


2.1 


1.9 


1.3 


Zti + nj 


5.6 


2.6 


2.6 


1.8 


2.3 


1.1 


1.1 


0.7 


Total 


185.3 


11.0 


10.5 


4.7 


87.2 


7.4 


4.4 


3.3 



Table 3: Same as Table §, but for the SSD + b channel. 



fb' 1 at 14(10) TeV. Therefore, we can conclude that the entire ranee from 500-1500 
GeV can be successfully probed at the LHC for = 0.2. 

• Stops decaying into higher neutralinos and charginos make the total rates distinctly 
better. This is governed by Higgsino couplings and is therefore most prominent for 
high tan/3 and low Aq. This effect is evident from the large event rates for point D. 
We can successfully probe this point in the SSD channel at 5<t with less than 1 fb^ 1 
data at both 10 and 14 TeV runs. 

• The trilepton final state occurs when the stop can decay into X21 X3 4 or 9- Therefore, 
points A and B show almost no signal and Point D has the largest signal in this 
channel. This advantage is largely lost for benchmark point E due to the kinematic 
suppression in the stop production process. 

• Reach for the LHC: Assuming the conservative case of (tan/3 = 5, Aq = 0), with 
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3/ 






14 TeV 








10 TeV 




Point 


Cut U 


Out 1 


Cut z 


Out 4 


Out (J 


Out 1 


Cut z 


Out 4 


A 


49.1 


2.8 


2.8 


0.0 


18.0 


1.1 


1.1 


0.0 


B 


2.0 


0.6 


0.6 


0.1 


1.1 


0.2 


0.2 


0.0 


U 


13.7 


y.i 


9.1 


1.1 


b.5 


3.7 


6.1 


O.o 




/IS 9 


9Q P. 

zy .o 


9Q n 
zy.u 


o.O 


9/1 1 
Z4b. 1 


i /i 

141.0 


141. 1 


41. D 


E 


9.3 


5.7 


5.5 


3.0 


3.4 


2.1 


2.1 


1.2 


tt + nj 


2.1 


0.0 


0.0 


0.0 


1.8 


0.0 


0.0 


0.0 


Wti + nj 


4.1 


2.5 


2.4 


1.0 


2.2 


1.4 


1.4 


0.6 


Zti + nj 


30.8 


20.7 


19.7 


2.7 


11.3 


7.3 


4.7 


1.1 


Total 


37.0 


23.2 


22.1 


3.7 


15.3 


8.7 


6.1 


1.7 



Table 4: Effect of cuts on signal and SM background cross sections (in fb) in the trilepton (31) 
channel at = 14, 10 TeV. Cut refers to all events passing the identification cuts. Cuts 1, 2 and 
4 are described in the text. Cut 4 is necessary to eliminate background from WZ + jets. 

10 fb" 1 luminosity, one can rule out \ e j^ greater than 0.007-0.045 (0.007-0.062) for 
stop masses between 500 and 1500 GeV at 95 % CL at <fs = 14 (10) TeV. A 5a 
discovery can be made in the same mass range for X' e ^ greater than 0.012-0.084 
(0.012-0.12). However, we observe that the reach in stop mass does not decrease 
monotonously with stop mass. The opening of new decay channels can improve 
detection considerably. The statements about minimum value of X' e ^ that can be 
probed are therefore dependent on the particular decays of the stop. We therefore 
tabulate the minimum values of X' e ^ for each benchmark point at 10 fb" 1 for both 
10 and 14 TeV in Table § 



Point 


14 TeV 


10 TeV 




95% CL 


3a 


5a 


95 % CL 


3a 


5a 


A 


0.007 


0.009 


0.012 


0.007 


0.009 


0.012 


B 


0.027 


0.037 


0.052 


0.029 


0.041 


0.059 


C 


0.026 


0.035 


0.048 


0.028 


0.038 


0.052 


D 


0.024 


0.032 


0.042 


0.027 


0.036 


0.047 


E 


0.045 


0.062 


0.084 


0.062 


0.087 


0.12 



Table 5: Values of minimum X' e ^ that can be ruled out at 95% CL, probed at 3er or ha with 10 
fb" 1 of data at yfs = 10, 14 for each of the benchmark points. The significance used is sj\fb where 
s is the signal and b is the background, s/b > 0.2 in all cases. 

In Figure ^, we present the effective mass distributions in SSD channel for all the 
benchmark points. Effective mass is denned as 

M eff = Yl ^ + zC ^ + P T ( 42 ) 

jets leptons 

The contributions from resonant stop production is superposed in the figures on the 
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SM backgrounds and also RPC superparticle production processes. The RPC contributions 
are much smaller and therefore do not provide a serious background to our signals. 
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Point C; 14 TeV 
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Figure 3: M e ff distributions at ^fs — 14 TeV. "SM" is the contribution to the background 
from Standard Model processes. "SUSY" refers to the contribution from R-conserving production 
processes. The inset in each figure contains the distribution for the signal alone. 



-13- 



4.2 Observability at the early run of 7 TeV 

The initial LHC run at y/s = 7 TeV will collect up to 1 fb _1 data. It will be difficult to 
observe RPV production of a 1 TeV stop at this energy. However, we can make useful 
comments for lower stop masses by looking at the SSD channel. We therefore look two 
benchmark points with low stop masses: the first is the 'Point A' described earlier and the 
second is similar to 'Point D' (with tan/3 = 10 and Aq = 0). Since the tt backgrounds are 
much smaller at 7 TeV, we relax the jet-pr cuts. The high scale parameters, stop mass at 
electroweak scale and cut-flow table for signal as well as background are given in Table ||. 

We conclude that we can rule out up to = 0.025 for a stop mass of 500 GeV 
at 7 TeV with 1 fb _1 data and a 5a discovery can be made at stop mass 500 GeV for 
^eff — 0.043. For the case tan /3 = 10 and Aq = 0, the lowest possible theoretically allowed 
stop mass (with m^j = 450 GeV) is 775 GeV and we can rule out up to X'^j = 0.054 with 
1 fb _1 data. 



Point 


(tan/3, A), ra ) 


M h 


Cut 


Cut 1 


Cut 2 


A 


(5,-1500,600) 


508 


283.3 


158.6 


147.5 


D' 


(10,0, 100) 


775 


70.0 


33.9 


32.0 


ti + nj 






116.0 


3.7 


3.5 


Wti + nj 






4.3 


2.3 


2.1 


Zti + nj 






1.8 


0.9 


0.8 


Total 




122.1 


6.9 


6.4 



Table 6: The benchmark points for studying RPV stop production and the effect of cuts on signal 
and SM background cross sections (in fb) in the SSD channel at y/s = 7 TeV. Cut refers to all 
events passing the identification cuts. All other cuts are described in the text, we do not apply Cut 
3. 



4.3 Differentiating from R-conserving signals 

We now address the question whether the signals we suggest can be faked by an R-parity 
conserving scenario in some other region(s) of the parameter space. One possible way that 
our signal may be mimicked is if a point in the mSUGRA parameter space (without RPV) 
gives similar kinematic distributions to any our benchmark points. More specifically, one 



may have a peak in the same region for the variable M e ff, defined in equation 4.2. 

For each of our benchmark points A-D, an M e ff peak in the same region requires 
the strongly interacting sparticles to have masses in the range already ruled out by the 
Tevatron data[34|. In particular, they require the gluino mass M g < 390 GeV. Thus the 
question of faking arises only for benchmark point E, which represents the highest mass 
where the signals rates are appreciable. 

We generate such a point (Point RC) with the parameters too = 300, nii/ 2 = 180, 
Aq = 0, tan/3 = 10, u > and the resultant sparticle masses for coloured particles are 
Mg = 465, Mq ~ 500 GeV. The M e jf distributions for point E and point RC is shown in 
Figure |]. We present the following results at 14 TeV as an illustration. Distributions at 
10 TeV are almost identical. 
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The missing Et distribution is also not a good discriminator under such circumstances, 
as can be seen from Figure |i[ This is because the neutrinos that contribute to missing- Et 
in the RPV case are highly boosted due to the large masses of the particles produced in 
the initial hard scattering. Thus, the Et spectrum is actually harder for the RPV case 
even though the RPC case has a stable massive LSP. However, as the resultant spectrum is 
quite light, the RPC production cross section (~ 40 pb) is about two orders of magnitude 
greater than the RPV case with A^ = 0.2 (~480 fb). Consequently, the rate of the SSD 
signals, for example, are much higher in the R-conserving scenario (~ 34 fb) as compared 
to those from point E (~ 14 fb). For values of A^ < 0.2, we can therefore make a reliable 
distinction simply based on the number of events expected in the SSD channel. 

Another possible discriminator is the charge asymmetry. In the SSD channel, one can 
look at the ratio of negative to positive SSD ^ ++ . The fraction of ds — > i* is more than 
the charge conjugate process ds — > t\ due to the difference in parton distributions of d 
and d in the proton. Therefore, one expects extra negative sign leptons than positive ones. 
Whereas in the RPC case, since most of the SSD contribution comes from gg production, 
we do not expect a large asymmetry. In our illustration, we see that this ratio is 2.7 (1.4) 
for the RPV (RPC) case. 




200 400 600 800 1000 1200 1400 1600 
M„ ff 




100 150 200 250 
Missing E T 



100 



Figure 4: The normalised effective mass (M e ff) and missing energy (fir) distributions in the 
SSD channel for Point E and an Point RC. Point RC has been generated using tuq = 300, = 
180; A = 0, tan/3 = 10 and u > 0. The gluino mass is 465 GeV. 



5. Non-x? LSPs 

For RPV models, the restriction of having an uncharged LSP no longer exits. A significant 
region of the mSUGRA parameter space with low uiq corresponds to a stau (f) LSP. With 
only A^-type couplings present, the stau can only decay via off-shell Xi an d t propagators 
into the four body decay (f — > rtds) if its mass, nif > mt op or via the five body decay 
(f — > rbWds) if mf < mt op where the top propagator is also off-shell. The four-body 



decays of the stau in lepton- number violating scenarios was calculated in [35]. Since the 
intermediate Xi i s °f Major ana character, we can always have one lepton of either sign 
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from LSP decay via the W from an on-shell or off-shell top. Thus, for various types of t 
decays, the following situations may arise: 

• For decays of ix^-type, we can still have same-sign dileptons with one lepton from 
top decay and the other from the decay of the LSP. 

• For decays of type b\f with \f ~^ W ± Xi + X, we have x? — > tt and the SSD come 
from W ± and LSP decay respectively. 

• For decays of the type t — > b\f with xf ~ * v tT + X, we may still get SSD from 
leptonic decay of the r in the f decay. If r-identification is used, final states of the 
type same-sign (r + e/fi) may be considered. 

• Since the stau has to decay via four- or five-body processes, it is possible that the 
lifetime of the f is large and it is stable over the length scale of the detector. In this 
case, it will leave a charged track like a muon and one can look at same-sign leptons 
with this "muon" as one of the leptons. It is also possible that the lifetime is large 
but the stau still decays within the detector. In this displaced vertex can be 
observed in the detector. 

We leave the detailed simulation of all scenarios of resonant stop production with stau LSP 
to a future study. 

Another possibility that arises with a large A 3 ' 12 is of having a stop LSP. In this case 
however, the decay will be almost entirely via the RPV ds di-jet channel. The over- 
whelmingly large dijet backgrounds at the LHC would most likely make this situation 
unobservable. 

6. Summary and Conclusion 

We have performed a detailed analysis of resonant stop production at the LHC, both 
for the 10 and 14 TeV runs, for values of the baryon number violating coupling A 3 ' 12 an 
order of magnitude below the current experimental limit. Benchmark points have been 
chosen for this purpose, which start just beyond the reach of the Tevatron and end close 
at the LHC search limit. We find that the same-sign dilepton final states, both with and 
without a tagged b, are most helpful in identifying the signal. The trilepton signals can 
also be sometimes useful, especially when decays of the resonant into higher neutralinos, 
the heavier chargino or the gluino open up. At 14(10) TeV, we can probe stop masses up to 
1500 GeV and values of A"jj down to 0.05(0.06) depending on the combination of various 
SUSY parameters. For cases of stop mass below a TeV, the effective mass distributions 
can enable us to distinguish between the resonant process and contributions from R-parity 
conserving SUSY processes. For higher stop masses, one has to rely on cross sections or 
the charge asymmetry. 

We have used a particular B-violating coupling, namely, X 312 - One can also have 
resonant stop production driven by A 3 ' 13 and A323. In either of these cases, one expects 
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a larger abundance of b quarks in the final state. However, there is a suppression in the 
production rates due to the b-distribution function in the proton. 

In conclusion, resonant stop production is a potentially interesting channel to look for 
SUSY in its baryon-number violating incarnation. Values of the B-violating coupling(s) 
more than an order below the current experimental limits can be definitely probed at the 
LHC, both at 10 and 14 TeV. If such interactions really exist, our suggested strategy can 
not only yield detectable event rates but also point towards resonant production as opposed 
to pair-production of SUSY particles. 
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